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Abstract 

This paper proposes a new mutual independence test for a large number of high di- 
mensional random vectors. The test statistic is based on the characteristic function of the 
empirical spectral distribution of the sample covariance matrix. The asymptotic distribu- 
tions of the test statistic under the null and local alternative hypotheses are established 
as dimensionality and the sample size of the data are comparable. We apply this test to 
examine multiple MA(1) and AR(1) models, panel data models with some spatial cross- 
sectional structures. In addition, in a flexible applied fashion, the proposed test can capture 
some dependent but uncorrelatcd structures, for example, nonlinear MA(1) models, multiple 
ARCH(l) models and vandermondc matrices. Simulation results are provided for detecting 
these dependent structures. An empirical study of dependence between closed stock prices 
of several companies from New York Stock Exchange (NYSE) demonstrates that the feature 
of cross-sectional dependence is popular in stock markets. 

Keywords: Independence test, cross-sectional dependence, empirical spectral distribu- 
tion, characteristic function, Marccnko-Pastur Law. 

1 Introduction 



A prominent feature of data collection nowadays is that the number of variables is comparable 
with the sample size. This is the opposite of the classical situations where many observations are 
made on low-dimensional data. Specific examples of such high dimensional data include microar- 
ray expression , images, multiuser dete c tion, climate data and financial data (see, for examp le, 



Donohol (j2nnd ). lFan. Lv and Qi 



(2011 



Johnstone (2001 



) . I Johnstone and Titteringtonl (|2009l )) 
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This type of data trend poses great challenges because traditional multivariate approaches do 
not necessarily work, which were establishe d for the ca s e of t he sample size n tending to infinity 



Anderson! ()1984l )). There hav e been a substantial 



and the dimension p remaining fixed (See 
set of r esearch work s dealing with high dimensional da t a (see , for example. 



Bai and Saranadasa 



(119961 ). iFan and Lvi (|2010l ). iHuang. Horowitz and Mai lj2008l ). iFan and Lil (|200lJ )). Measuring 
mutual dependence is important in time series analysis and cross-sectional panel data an alysis. 
While s erial depend ence can be characterized by the general spectral density function (see iHong 



l|l998l ): 



Hon 



3 ( 19991 



)), mutual dependence is difficult to be described by a single criteria. This 



paper proposes a new statistic to test mutual dependence for a large number of high dimensional 
random vectors, including multiple time series and cross-sectional panel data. 

Suppose that {Xji,j = l,...,n;i = l,...,p} are real-valued random variables. For 1 < 
i < p, let Xj = (Xu, • • • , X n i) T denote the i— th time series and xi, • • • , x p be a panel of p time 
series, where n usually denotes the sample size in each of the time series data. In both theory 
and practice, it is not uncommon to assume that each of the time series (Xu,X2i, ■ ■ ■ ,X n i) is 
statistically independent, but it may be unrealistic to assume that xi, X2, • • • , x p are independent 
or even uncorrelated. This is because there is no natural ordering for cross-sectional indices. 
There are such cases in various disciplines. In economics and finance, for example, it is not 
unreasonable to expect that there is significant evidence of cross-sectional dependence in output 
innovations across p countries and regions in the world. In the field of climatology, there is also 
some evidence to show that climatic variables in different stations may be cross-sectionally 
dependent and the level of cross-sectional dependence may be determined by some kind of 
physical distance. Moreover, one would expect that climatic variables, such as temperature and 
rainfall variables, in a station in Australia have higher-level dependence with the same type of 
climatic variables in a station in New Zealand than those in the United States. 

In such situations, it may be necessary to test whether xi,X2, • • • , x p are independent before 
a statistical model is used to model such data. In the econometrics and statistics literature, sev- 
eral papers have basically considered testing for cross-sectional unc orrelatedness fo r the residuals 
involved in some specific regression models. Such studie s include Pesaranl (2004) for the para- 



metric linear model case. 



Hsiao. Pesaran and Pick! (120071 ) for the parametric nonlinear case, and 



Chen. Gao and Lil (20091 ) for the nonparametric nonlinear case. Other related papers include 



Su and Ullah 



2009j) for testing conditional uncorrelatedness through examining a covariance 



matrix in the case where p is fixed. As the main motivation of this paper, we will propose using 
an empirical spectral distribution function based test statistic for cross-sectional independence 
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of Xi,X 2 , • • • ,x p . 
The aim is to test 



H : xi 

where Xj 



, x p are independent; against Hi : xi, 



x p are not independent, 



(1.1) 



(Xu, . . .,X ni ) T for i = 1, . . . ,p. 

In time series analysis, mutual independence test for multiple time series has long been 
of interest. Moreover, time series always display various kinds of dependence. For example, 
an autoregressive conditional heteroscedastic (ARCH(l)) model involves a martingale difference 
sequence (MDS); a nonlinear moving average (MA) model is not a MDS, but its autocorrelations 
are zero; a linear moving average (MA) model and an autoregressive (AR) model are both 
models with correlated structures. In this paper, we also employ the proposed statistic to test 
dependence for multiple time series. 

Section 8.5 of I Anderson! ()1984l ) also consider a similar problem but with fixed dimensions. 
His problem and approach are as follows. Let the pm-component vector x be distributed ac- 
cording to N(fi, E). Partition x into m subvectors with p components respectively, that is, 
(h^, • • • , ) T . The vector of means /j, and the covariance matrix S can be partitioned respec- 



tively, (a*iV" ,/Xp) T , and 



^21 ^22 



(1.2) 



/ 



The null hypothesis to be tested is that the subvectors hi, I12, . . . , h p are mutually indepen- 
dently distributed. If hi, Ii2, . . . , h p are independent subvectors, 



E(xi - fii){xj - fij) 



0, Vi ^ j. 



(1.3) 



Thus, the null hypothesis is equivalent to testing Ho : Sj,- = 0, Vz 7^ j. This can be stated 
alternatively as the hypothesis that 5] is of the form 



En 
£22 










(1.4) 



3 



To tackle the problem above, draw n observations from the population x to form a sample 
covariance matrix. The likelihood ratio criterion proposed is 

|Q|3» 



A 



nf=iicy*' 



(1.5) 



where 



Q 



Qn Q12 
Q21 Q22 



Q,2p 



(1.6) 



i Qpi Qp2 • • • Qpp 1 

in which Qjfc is the sample covariance matrix of the random vector hj and h^. 

Our approach essentially uses the characteristic function of the empirical spectral distri 
bution of sample covariance matrices in large random matrix theory. When x^ , • • • , Xp ar< 
mutual independent, the limiting spectral distribution (LSD) of the cor responding sample co- 
varian c e matrix is the Marcenko -Pastur (M-P) law (see, for example, 



Marcenko and Pastur 



(119671 ) . iBai and Silversteinl (12009)1. From t 



evidence of dependence. Indeed, 



l is poi nt, a ny deviation of the LS D from M-P law is 



Silversteinl (j 19951 ) and lBai and Zhoul (|2008l ) report the LSD of 



the sample covariance matrix with correlations in columns and it is different from the M-P law. 
In dependent but uncorrelated cases, the LSD may also be the M-P law, e.g. ARCH(l) model. 
However, the proposed test can still capture the dependence by utilizing the nonzero correlation 
between high-order series. Unlike the Anderson's test, we need not re-draw observations from 
the set of vectors of xi , • • ■ , x p due to the high dimensionality. 

The rest of the paper is organized as follows. Section 2 introduces the proposed test statistic 
and some related large dimensional random matrix results. We present the asymptotic distri- 
butions of the test statistic under both the null and local alternative hypotheses in Section 3. 
Moreover, its applications to some panel data models with spatial dependent structures are also 
illustrated. Section 4 studies a general panel data model whose dependent structure is different 
from that investigated in Section 3 and develops the asymptotic distribution of the proposed 
statistic under the null hypothesis for this model. Several simulated results which demonstrate 
the effectiveness of the proposed test under the two dependent structures are shown in Section 
5. Furthermore, this section illustrates simulation results for some dependent but uncorrelated 
structures, though we have not developed asymptotic theory for these kinds of models. An 
empirical analysis of daily closed stock prices from NYSE is provided in Section 6. We conclude 
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with a discussion in Section 7. All mathematical proofs are relegated to an appendix available 
from a supplementary document. The full version of this paper is available as a working paper 
at http://www.jitigao.com 



2 Theory and Methodology 

Before we establish the main theory and methodology, we first make the following assumptions: 

Assumption 1. For each i = 1, . . . ,p, Yu, ■ ■ ■ ,Y n i are independent and identically distributed 
(i.i.d) random variables with mean zero, variance one and finite fourth moment. When Yji 's are 
complex random variables we require EY^ = 0. Let Xj = T^Vj with y« = (Yu, ■ ■ ■ ,Y n i) T and 

1/2 

T n being a Hermitian square root of the nonnegative definite Hermitian matrix T„ . 

Assumption 2. p = p(n) with p/n — > c 6 (0, oo). 

We stack p time series Xj one by one to form a data matrix X = (xi, • • • ,x p ). Moreover 
denote the sample covariance matrix by 

A n = ix*X 

n 

where X* is the Hermitian transform of the matrix X. The empirical spectral distribution (ESD) 

of the sample covariance matrix A n is defined as 

v 



F A ^x) = -^I(X l <x), (2.1) 

where Ai < A2 < . . . < A p ar6 6ig6nvaliiGS of A n . 

It is well-known that if independent and p/n — > c £ (0, 00) then F An (x) 

con yerges with probabil i ty one to the Marcenko-Pastur Law F c (x) ( see Marcenko-Pastur (1967) 



and 



Bai and Silversteml (|2009l )) whose density has an explicit expression 



fc(x) 



0, otherwise; 

and a point mass 1 — 1/c at the origin if c > 1, where a = (1 — yfc) 2 and b = (1 + yfc) 2 . 

When there is some correlation among xi, • • • , x p , denote by T p the covariance matrix of the 
first row, yj , of X. Then, under Assumption 1, when F Tp (x) —> H(x), F An converges with 
pro bability one to a n on random distribution function F C,H whose Stieltjes transform satisfies 



(see ISilversteinl ()1995l )) 

m{z) = [ — 1 — f - dH(x), (2.3) 

J x(l — c — czm(z)) — z 
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where the Stieltjes transform mg for any c.d.f G is defined as 

m G ( z ) = J —L-dG(X), hnz > (2.4) 

and G can be recovered by the inversion formula 

1 f X2 

G{[xi, X2]} = — lim / Zm{mG{x + ie))dx, 
vr e^0+ J X1 

where x\ and x% are continuity points of G. 

Moreover, equation (|2.3[) takes a simpler form when F C,H is replaced by 

F c ^ = (l-c)/ [0>oo] + C F c ^, (2.5) 

which is the limiting ESD of A n = ^XX*. Its Stieltjes transform 

1 - c 

m{z) = h cm(z) (2-6) 

has an inverse 

z = z (m) = — + c[ — - — dH(x). (2.7) 
m J 1 + xm 

The construction of our test statistic relies on the following observation: the limit of the ESD 
of the sample covariance matrix A n is the M-P law by (I2.2j) when xi , • • • , x p are independent 
and satisfy Assumption 1, while the limit of the ESD is determined from (|2.3|) when there is 
some correlation among xi , • • • , x p with the covariance matrix T p different from the identity 
matrix. Moreover, preliminary investigations indicate that when xi,-- - ,x p are only uncorre- 
lated (without any furth er assumptions) , the limit of the ESD of A n is not the M-P law (see 



Ryan and D ebb an 



( 20091 )). These therefore motivate us to employ the ESD of A n , F n (x), as 
a t est statistic. There is no central limit theorem for (F An (x) — F c,H (x)), however, as argued 
by lBai and Silversteinl ( 20041 ) . We instead consider the characteristic function of F An (x). 



The characteristic function of F n (x) is 



f 1 ' 

Sn (t) 4 J e itx dF An (x) = -JZe itx ', (2i 



i=l 

where Ai, i = 1, . . . , p are eigenvalues of the sample covariance matrix A n . 
Our test statistic is then proposed as follows: 

fT 2 



M n = I \s n (t)-s(t)\ 2 dU(t), (2.9) 

Ti 
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where s(t) := s(t,c n ) is the characteristic function of F Cn (x), obtained from the M-P law F c (x) 
with c replaced by c n = p/n, and U (t) is a distributional function with its support on a compact 
interval, say [Ti,T2]. 

To develop the asymptotic distribution of M n under a local alternative, the following as- 
sumption is needed. 

Assumption 3. Let T p be apxp random Hermitian nonnegative definite matrix with a bounded 
spectral norm. Let yj = wjTp , where T^/ 2 is the p x p Hermitian matrix that satisfies 
(T)!/ 2 ) 2 = T p and wj = (Wji, • • • , Wj p ),j = 1, . . . , n are i.i.d random vectors, in which Wji,j < 
n,i <p are i.i.d with mean zero, variance one and finite fourth moment. 

The empirical spectral distribution F Tp of T p converges weakly to a distribution H on [0,oo) 
as n — > oo; all the diagonal elements of the matrix T p are equal to 1. 

Note that under Assumption [3l A n becomes Tp^ 2 W*WTj!/ 2 , where W = (wi, • • • ,w n ) T . 
The assumption that all the diagonal elements of T p are equal to 1 is used to guarantee that 
EXji = 1. Roughly speaking, this constraint is to ensure that the data are weakly stationary. 
However, we would like to remark that such a constraint could be removed if we knew that T p 
was not a diagonal matrix. Under Assumption |3J when T p = I p , the random vectors xi, . . . , x p 
are independent and when T p ^ I p , they are not independent. For convenience, we name this 
dependent structure as 'linear dependent structure'. 

To develop the asymptotic distribution of the test statistic, we introduce 

G n (x) = p[F An {x) - F Cn {x)}. (2.10) 

Then, p(s n (t) — s(t)) can be decomposed as a sum of the random part and the non-random 
part as follows: 

p(s n (t) - s(t)) = J e itx dG n (x) 

= J e itx d(p [F An {x) - F Cn ' Hn {x)}) + J e itx d(p [F°»> H »(x) - F Cn {x)]) , (2.11) 
where F Cn,Hn is obtained from F C,H with c and H replaced by c n = p/n and H n = F Tp . 

Theorem 1. Let Assumptions^ and® hold. 

1) Suppose that Assumptions^ and the following conditions are satisfied: 

n 

- £ e* Ty\m{ Zl )T p + i)" 1 T^e* (m(z 2 ) T p + i)" 1 -> h 1 (z 1 ,z 2 ) (2.12) 
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and 

n 

- £ e* Ty 2 {m{z) T p + i)" 1 T^e.e* T^irn^z) T p + i)~ 2 T^e, -> /i 2 (z), (2.13) 



n 
i=l 



where e* zs i/ie n- dimensional row vector with the i-th element being 1 and others 0. 

Then, the scaled proposed test statistic p 2 M n converges in distribution to a random variable 
R of the form: 

R= [ \\V(t) + 5i(t)| 2 + |Z(t) + 5 2 (t)\ 2 )dU(t), (2.14) 
where (V(t),Z(t)) is a vector of Gaussian processes and Sj(t), j = 1,2 are defined as 

<5i(£)=lim cos(tx)dp(F c "' Hn (x) - F c "(x)), (2.15) 

n— >oo J ' 

fatt) = lim sin(tx)dp(F c "' Hn (x) - F c "(x)). (2.16) 

ra— >oo J ' 

The mean and variance of (V(t),Z(t)) are specified as follows. If Xn is real, then 

1 f cjm 3 (z)r 2 (l + rm(z))- 3 dff(r) 

^ (t) = 2^ (1 - c/ m^)r 2 (l + m(,))-W(r))^ Z 

ggn - 3 i cos(M cnfWWz) dz 

Replacing cos{tz) in EV(t) by sin(tjz) yields the expression of EZ(t). Also 

1 / / cos(tjZi)sin(tf l Z2) d d 



^ (*r, \ r,, \\ 1 7 7 cos{tjZi)sin(thZ2) a . . ct , . , 

Cm; F =--^i i Jl m( Z2 ) m (z 1 )dz 1 dz 2 



c(EX^ — 3) I I d? 

■ cos{tjZi)sin(t h z 2 )- — ^[m(zi)m(z2)hi(zi, z 2 )\dz\dz 2 . 



4vr 2 /yj ^ dz\dz 2 

If X\\ is complex with EX^ = 0, then 

. . £|X n | 4 -2 / . cm 3 (z)h 2 (z) 

w 2vri J 7 v 7 1 - cjm 2 (z)T 2 dH(T)/(l + Tm(z)) 2 

and the covariance 

n fvu \ vs. \\ 1 I I cosCtjZ^sinCthZz) d d . . 

V / 47r z y 7l y 72 (m(zi) - m{z 2 )Y dz 2 dz\ 

c{E\X\\fi — 2) / / . . . . cl 2 

1 cos[tjZ\)sin{t h z 2 )- — - — \jn(zi)m(z 2 )hi(zi, z 2 )\dz 1 dz 2 . 



4vr z J 7l J 72 dz\dz 2 

The covariances Cov(Vj,Vh) and Cov(Zj, Zh) are similar to Cov(Vj, Zh) except replacing 
(cos(tjz), sin(th,z)) by (cos(tjz),cos(thz)) and (sin(tjz), sin(thz)) respectively. 
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The contours 7, 71 and 72 above are all closed and are taken in the positive direction in the 
complex plane, each enclosing the support of F C,H . Also 71 and 72 are disjoint. 

2) Under the null hypothesis Ho, the scaled statistic p 2 M n then converges in distribution to 

oT 2 



Ti 



(\v(tr + \z( t )\')du(t), 



(2.17) 



where the distribution of (V(t),Z(t)) can be obtained from (V(t), Z(t)) with H{t) becoming 
the degenerate distribution at the point 1, m(z) being the Stieltjes transform of the M-P law, 

M*l> Z ?) = (m(*i)+l)W*2)+l) a7ld = (m(4+l) a • 



Remark 1. Assumption^ assumes that all the entries o/xj are identically distributed. It is of 
practical interest to consider removing the identical distribution condition. Instead of assuming 
identically distributed entries for Xj, we need only to impose the following additional assumptions: 
for any k = 1, . . . ,p; j = 1, . . . ,n, EX jk 



0, EX% 



1, sup EXj k < 00 and for any 7] > 0, 



* n P 



n A np 



0. 



(2.18) 



j=i k=i 



A caref ul checking on t h e argu ments of Theorem 1.1 o nBai and Silverstein 1200 4 ) and The- 
orem 1.4 of \Pan and Zhou $200a ) indicates that Lemma 6 (of Appendix A of the supplementary 
material) still holds and hence TheoremUl holds under $2.18\) . For the expressions of the mean 
and covariance of the asymptotic random vector, we substitute the fourth moment E\Wii\ 4 in 



Lemma 6 with the average of all the fourth moments of all the entries, i.e. 



np 



Conditions (|2.12p and (|2.13p can be removed if E \W\i\ = 3 in the real-number case or 
E \ W\i\ = 2 in the complex-number case (see, for example, iBai and Silverstein! ( 2004! ) ) . The 



second part of the above theorem is concerned with asymptotic distributions of the test statistic 

under a local alternative hypothesis, i.e., Assumption [3j With respect to Assumption [3l we 

would make the following comments, which are useful in the subsequent application section. 

If yj = wJC, where C is any q x p nonrandom matrix and Wj, j = 1, . . . , n are i.i.d. gx 1 

random vectors with their respective entries being i.i.d random variables, then Theorem [1] is 

still applicable. This is because ^XX* in this case becomes ^WCC*W* and the nonnegative 

1/2 1/2 1/2 

definitive matrix CC* can be decomposed into CC* = T q T q , where T q is a q x q Her- 

mitian matrix and W = (#1, . . . , w n ). Note that the eigenvalues of ^XX* differ from those of 
I T i/2 w * WT i/2 by | p _ ? | zerog _ Thug) we may ingtead resort tQ CLX of lTq /2 W*WTg /2 . 
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3 Applications in multiple MA(1), AR(1) and spatial cross- 
sectional dependence structures 



This subsection is to explore some applications of the proposed test. 

Example 3.1. Consider a multiple moving average model of order 1(MA(1)) of the form: 

v t = z t + ipzt-i,t = l,...,p, (3.1) 

where l^l < 1; = (Zu, . . . ,Z n t) T is an n-dimensional random vector of i.i.d. elements, each 
of which has zero mean and unit variance; and vt = (Vu, ■ ■ . , V n t) T . Denote by Vj and Zj 
respectively the j-th rows of V = (V jt )nxp and Z = (Z jt ) nx ( p+1 y 
For each j = 1, . . . , n, the MA(1) model (13. ip can be written as 



j 



zjc, 



(3.2) 



where 



t 





• 


■ 





1 




• 


• 








1 


i> ■ 


• 











l v 

1 



(3.3) 



\ / (p+l)xp 

From Assumption [3l the last paragraph of the preceding subsection and Theorem [H our test 
is able to capture the dependence of vi, • • ■ , v p as n and p go to infinity in the same order. 

Example 3.2. Consider a multiple autoregressive model of order 1(AR(1)) of the form: 

1 



v 4 = 0v t _i + z t , t = l,...,p, v 



=zo, 



(3.4) 



where \4>\ < 1; for any t = 0, 1, . . . ,p, Zt = (Zu, ■ ■ ■ , Z n t) T is an n-dimensional random vector 
with i.i.d. elements, each of which has zero mean and unit variance; and v t = (Vu, ■ ■ ■ , V n t) T • 
Denote the j-th rows of V = (Vjt) nx ( P +i) and Z = (^jt) n x(p+i) as Vj and Zj respectively. 
For each j = 1, . . . , n, the AR(1) model (|3.4|) can be written as 



V 



ZjD, 



(3.5) 
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where 



D 



1 -0 (-0 

o l -4 

1 



( _0)p-2 ( _^ r i/yr^2 

P-3 (_0)p-2/yf302 



V 



o o 










1 





1A/1 



(3.6) 



By Theorem [IJ we can also apply the proposed test M n to this AR(1) model as well. 

Example 3.3. We now consider a panel data case. Let {vji : i = 1, . . . ,p;j = 1, . . . , n} be the 
error components in a panel data model. They may be cross-sectionally correlated. In panel 
data analysis, it is of interest to consider the cross-sectional independence hypothesis, i.e. 

Hoo : Cov(vji,Vjh) = for all j = 1, . . . , n and all i ^ h; 

against 

Hn : Cov(vji,Vjh) 7^ for some j and some i ^ h. (3-7) 

Under the assumption that {vji : i = 1, . . . ,p;j = 1, . . . , n} are normal distributed, this hypoth- 
esis is equivalent to the independence hypothesis that 

Ho : Vi, . . . , v p are independent; against Hi : vi, . . . , v p are not independent, (3.8) 

where Vj = (v u , . . . ,v ni ) T , i = l,...,p. 

Modern panel data literature has mainly adopted two different approaches to model error 
cross-sectional dependence: the spatial approach and the factor-structure approach. For the 
spatial approach, there are three popular spatial models: Spatial Moving Average (SMA), Spa- 
tial Auto-Regressive (SAR) and Spatial Error Components (SEC) processes. They are defined 
respectively as follows: 

p 

SMA : Vji = ^UikSjk + £jh ( 3 -9) 



k=l 

p 



SAR : Vji = '^uJikVjk + £ji, 
k=i 
p 

SEC : Vji = uJikCjk + 



(3.10) 
(3.11) 



fc=i 
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where Uik is the i-specific spatial weight attached to individual k; {eji : i = 1, . . . ,p;j = 1, . . . , n} 
and : i = 1, . . . , p; j = 1, . . . , n} are two sets of i.i.d. random components with zero mean and 
unit variance, and : i = 1, . . . ,p; j = 1, . . . ,n} are uncorrelated with {eji,i = 1, . . . ,p;j = 
l,...,n}. 

Denote the j-th row of V = (vji) nXp , e = (eji) nxp and £ = (^i)nxp by vj, ej and £j 
respectively. Set u> = {cJik)pxp- Model SMA (j3.9|) may be rewritten as vj = eJ(cj T + I p ), Vj = 
1, . . . , n and hence T p = (u+I p )(u T +I P ). For model SAR (fHTTUD . assume that u— I p is invertible. 
We then write vj = eJ(uj t — Vj = 1, . . . ,n. Hence T p = (a; — I p ) _1 (a; T — I p ) _1 . Therefore 
the test statistic M n can be used to identify whether vi , • • ■ , v p of models (|3.9p and (|3.10p are 
independent. Hence it can capture the cross-sectional dependence for the SMA model and SAR 
model 

As for the SEC model defined in ()3.11j) . whether the statistic M n can detect the dependence 
structure of the SEC model relies on the properties of a sample covariance matrix of the form: 

1 
— 

n 

where $ = . . . ,£ P ) T and e = ( ej , . . . , e p ) T . 



B n = -(u£ + e)(uZ + e) T , (3.12) 



Under the null hypothesis Ho, iDozier and Silverstenl (|2007l ) provides the LSD of the matrix 
B n whose Stieljes transform is 

™(*) = I —x ,? H{ - X ), ^^ ; , (3-13) 

J i ■ - ( \ - (1 + crh(z))z + 1 - c y ' 

J l+crn(z) v \ J/ • 

IT 

where H(x) is the limit of F~ ££ . 

With this result, we know that the LSD of the matrix B n is not the M-P law. In view of 
this, the proposed test M n should capture the dependence of the SEC model (|3. 1 1 j) in theory. 
However, the asymptotic distribution of the proposed test statistic M n for this case will need to 
be developed in future work. 



4 A general panel data model 

Note that the proposed test is based on the idea that the limits of ESDs under the null and local 
alternative hypotheses are different. Yet, it may be the case where there exists some dependence 
among the set of vectors of vi, • • • , v p but the limit of the ESD associated with such vectors is 
the M-P law. Then a natural question is whether the statistic M n works in this case. We below 
investigate a panel data model as an example. 
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Consider a panel data model of the form 

Vij = Eij + -^ u ii « = !)••• >P; j = 1, . . . , n, (4.1) 

where {£ij,i = 1, . . . ,p; j = 1, . . . , n} is a sequence of i.i.d. real random variables with Ee\\ = 
and Ee\ x = 1, and {iti,i = 1, . . . ,p} are real random variables, and independent of {£ij,i = 
1, = l,...,n}. 

For any i = 1, . . . , p, set 

Vj = (v a , . . . ,v in ) T . (4.2) 

The aim of this section is to test the null hypothesis specified in ()3.8|) for model ( ]4.ip . 
Model (|4.ip can be written as 

v = e + ue T , (4.3) 

where v = (vi, . . . , v p ) T , u = (-^m 15 . . . , -^=u p ) T and e is p x 1 vector with all elements being 
one. 

Consider the sample covariance matrix 

S n = ivv T = -(e + ue T )(e + ue T ) T . (4.4) 
n n 

By Lemma [5] in the appendix and the fact that rank(ue T ) < 1, it can be concluded that 
the limit of the ESD of the matrix S is the same as that of the matrix j^e 1 \ i.e. the M-P 
law. Even so, we still would like to use the proposed statistic M n to test the null hypothesis of 
mutual independence. However, this model does not necessarily satisfy Assumption [1] because 
the elements of each vector Vj are not independent and they include the common random factor 
Ui. As a consequence, Theorem Q] thus can not be directly applied to this model. Therefore, we 
need to develop a new asymptotic theory for the proposed statistic M n for this model. 

Theorem 2. For model |^.1[ ), in addition to Assumptions^ and\^ we assume that 

p 



i?||u|| 4 < oo and \e 



^2( u l - u ){ u ) - «) 



as n — >■ oo, (4-5) 



where u is a positive constant. 

Then, the proposed test statistic p 2 M n converges in distribution to the random variable R2 
given by 

R 2 = t (\W(t)\ 2 + \Q(t)\ 2 )dU(t), (4.6) 
Jti 
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where (W(t),Q(t)) is a Gaussian vector whose mean and covariance are specified as follows: 

^ rTT r/M c / , , m 3 (z)(l + m(z)) 

E\W(t)} = & cosftjz) ,. — \\\ ~ v ,^ dz 

L WJ 2vriJ 7 3 ; ((l + m(z)) 2 -cm 2 (z)) 2 

c(EY& - 3) / , m 3 (z) 

2iri % COS{tjZ) (1 + m{z)Y - cw?{z) dZ 



and 



c f . . m(z) 

H COS(tjZ)—- : — TTo n . . , az 

2vriJ 7 ° z((l + m(z)) 2 -cm 2 (z)) 



1 jf / cos(tjZi)sin(thZ2) d ^ d 



Cov(W(tj),Q{t h )) = -—r f d> \ m(zi)—m(z 2 )dz 1 dz 2 

2vH / 71 / 72 -rn(z 2 )) 2 dz x dz 2 

c(EYA-3) f f , , , (i 2 r m(z 1 )m(z 2 ) , 

7^ ? f cos{tjz{)cos(t h zi) [j— — p rw- - — p-^ dzid2&. 

4vr 2 / T1 / 72 dzidz 2 (1 + + m(^2)) 



Replacing cos(tjz) in E[W(tj)] by sin(tjz) yields the expression of E [Q(tj)]. The covariances 
Cov(W(tj),W(th)) and Cov(Q(tj),Q(th)) are similar except replacing sin(thz) and cos(tjz) by 
cos(thz) and sin(tjz) respectively. The contours in f^, 7j ) and both enclose the interval 

[(1 — y/c) 2 + 2cu, (1 + \/c) 2 + 2cu]. Moreover, the contours 71 and 72 are disjoint. 

Remark 2. When u±, ■ ■ ■ ,u p are independent and hence v±, ■ ■ ■ , v p are independent, condition 
\4-5\) is true. 

In view of Theorem [21 we see that the proposed test statistic M n still works mainly due to 
the involvement of the last term on the right hand of (|2.1ip . 



5 Small sample simulation studies 



This section provides some simulated examples to show the finite sample performance of the 



proposed test. In addition, we also 



a likelihood ratio test proposed by 



compare t 



le pe rformance of the proposed test with that of 



Andersonl (119841 ) . Simulations are used to compute and find 



the empirical sizes and powers of the proposed test and hence evaluate the performance of the 
test. To show the efficiency of our test, the two kinds of dependence structures investigated in 
Section 3 and 5 are detected, such as multiple MA(1) and AR(1) model, SMA and the general 
panel data model. 
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5.1 Empirical sizes and power values 

First we introduce the method of calculating empirical sizes and empirical powers. Let z\„ and 
z 1 _i a be the 100(|a)% and 100(1 — quantiles of the asymptotic null distribution of the 

test statistic M n respectively. With K replications of the data set simulated under the null 
hypothesis, we calculate the empirical size as 

U of M» > Zl _i a or M« < zi a } 
a = ^ 2 —, (5.1) 

where represents the values of the test statistic M n based on the data simulated under the 
null hypothesis. 

In our simulation, we choose K = 1000 as the number of repeated simulations. The sig- 
nificance level is a = 0.05. Since the asymptotic null distribution of the test statistic is not a 
classical distribution, we need to estimate the quantiles zi a and z 1 _i a . Naturally, we do as 
follows: generate K replications of the asymptotic distributed random variable and then select 
the (K^a)-th smallest value zi a and (K^a)-th largest value z 1 _i a as the estimated 100(^a)% 
and 100(1 — \ol)% quantiles of the asymptotic distributed random variable. 

With the estimated critical points Zi a and z 1 _i a under the null hypothesis, the empirical 
power is calculated as 

A {0 of > %_i Q or < zij 

P= 1 2 ~, (5-2) 

where represents the values of the test statistic M n based on the data simulated under the 
alternative hypothesis. 

5.2 Comparisons with the classical likelihood ratio test 

For the proposed independence test, we generate n numbers of p-dimensional independent and 
identical distributed random vectors {yj}™ =1 , each with the mean vector p and the covariance 
matrix S. Under the null hypothesis, {yj}" = i are generated in two scenarios: 

1. Each Wj is a p-dimensional normal random vector with the mean vector P and the co- 
variance matrix S = I p for j = 1, . . . n, y,- = T p Wj with T p = l p ; 

2. Each Wj consists of i.i.d. random variables with standardized Gamma(4,2) distribution, 
so they have zero means and unit variances; for j = 1, . . . n, Vj = T p Wj with T p = l p . 

Under the alternative hypothesis, we consider the case: T p ' = (y/oMlp, V0M1 P ), where 
l p is a p-dimensional vector with 1 as entries. In this case, the population covariance matrix of 
Vj is £ = 0.95I P + 0.051 p l p , which is called the compound symmetric covariance matrix. 
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For normal distributed data, the fourth moment of each element is E'lXnj 4 = 3; for stan- 
dardized Gamma(4,2) distributional data, -ElXiiJ 4 = 4.5. This can be calculated by a formula 
as follows: the A;-th moment of X\\ which is Gamma(a, A) is 

E\X u \ k = (a + fc - 1 V;" (a + 1)a . (5.3) 
A fc 



Anderson! (|1984l ) provides a likelihood ratio criterion (LRT) to test independence for a fixed 



number of fixed dimensional normal distributed random vectors. We compare it with the pro- 
posed test. 

Under the null hypothesis, the distribution of L is the distribution of L2-L3 • • • L p , where 
L2, . . . ,L p are independently distributed with having the distribution of U m ^k-i)m,n-i~(k-i)m- 
Furthermore, for any k = 2, . . . ,p, as n — > 00, — (n — | — ^)log [ y m ./fc-iw 1 p-i-(fc- i)m] has x 2 
distribution with (k — l)m 2 degrees of freedom (see Section 8.5 of lAndersonl (|1984l )). 

From the construction of the LRT test, we can see that the LRT utilizes additional n obser- 
vations of the random vectors xi, . . . ,x p under investigation, while the proposed test does not 
need this information. However, we can choose m = 1 and apply LRT for testing the indepen- 
dence of the random vectors Xi, . . . ,x p , where for any i = 1, . . . ,p, the elements of the vector 
Xj consist of its n observations. Hence the LRT test can test independence for the number of p 
random vectors with dimension n by choosing m = 1. 

Tables [T] and [5] show the empirical sizes and power values of our proposed test and the LRT 
test for normal distributed random vectors respectively. From Tables Q] and [21 we can see that 
the LRT test does not work when p and n are both large while the proposed test possesses good 
performance when p and n go to infinity at the same order. The LRT test is only applicable to 
the case where p is fixed and n tends to infinity. From Table [2 it can be seen that the LRT 
fails when p is comparable with n. When the difference between p and n are large, the sizes 
and the powers of the proposed test become worse. This is because our test is proposed under 
the case that p and n are in the same order when they approximate to infinity. The proposed 
test also works well for gamma random vectors while the LRT test is not applicable to gamma 
case, since, in theory, LRT test is provided for normal random vectors. Table [3] provides the 
empirical sizes and empirical powers of the proposed test for the gamma case. In our simulation, 
we choose p, n = 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100 for the proposed test and the LRT test. 
The significant level a is chosen as 0.05. In each case, we run K = 1000 repeated simulations. 
Reports for empirical powers show that the proposed test can check independence for both 
normal and gamma vectors well. Moreover, the empirical powers converge to 1 as n,p-> 00. 
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5.3 Multiple MA(1), AR(1) and SMA model 

Consider multiple MA(1) model 



Vt = z t + ipz t -i ) t=l,...,p. (5.4) 

We choose tp = 0.5 and the simulation results in Table H] show that the proposed test performs 
well for this model. 

Consider multiple AR(1) 

v t = <pwt-i + Zf, t = 1, . . . ,p, v = -, 1 z . (5.5) 

V 1 - 

Let ^> = 0.5. The empirical powers for this model are provided in Table As n and p increase 
in the same order, the empirical power tends to 1. 

As for the Spatial Moving Average (SMA) model, i.e. 

p 

Vji = ^2 ^ikZjk + (5.6) 

k=l 

generate Ejj. *~ normal(l, 1), Vj = 1, . . . , n; k = 1, . . . ,p. 

Applying the proposed statistic M n for the sample matrix - V* V, the empirical power values 
given in Table [6] show that M n performs well for capturing the cross-sectional dependence for 
SMA model. 



5.4 The general panel data model 

We examine the finite sample performance of the proposed test for the general panel data model 
(BID, i.e. 

Vij = Eij + —/= u ii i = 1) • • • >2>; J = 1, . . . , rx, (5.7) 

where {£ij,i = 1, . . . ,p; j = 1, . . . , n} is a sequence of i.i.d. random variables with Ee\\ = and 
Ee\y = 1, and {ui,i = 1, . . . ,p} are independent of = 1, . . . ,p; j = 1, . . . , n}. 

Under the null hypothesis, we generate Ui l ~ d normal ( 1, 1), i = l,...,p and under the 
alternative hypothesis, we generate the data from u = (-j=m, -^ u 2, • • • , ~jf; u p) ~ ^p^^p } 
where S = TT T and T is a p x p matrix with its elements being generated by * U(0, 1), 

= 1,... ,p. 

The simulation results for the empirical sizes and power values given in Table [7] show that 
the proposed test can capture the dependence for the general panel data model (I4.ip . 
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5.5 Some other time series models and Vandermonde matrix 

Dependent structures of a set of random vectors are normally described by non-zero correlations 
among them, such as the linear dependent structure developed in Section 3. However, there are 
some data that are not independent but uncorrelated. We consider three such examples and test 
their dependence by the proposed test although we have not developed an asymptotic theory 
for each of such local alternatives. 

5.5.1 Nonlinear MA model 

Consider a nonlinear MA model of the form 

R tj = Z t -ijZ t -2,j(Z t -2,j + Z tj + 1), t = l,...,p; j = 1, . . . , n; (5.8) 

where zt = (Zn, . . . , Zt n ) is an n-dimensional random vector with i.i.d. elements, each of which 
has zero mean and unit variance, and r^ = (Rti, • • • , Rtn)- 

For any j = 1, . . . , n, the correlation matrix of (Rji,Rj2, ■ ■ ■ , Rj P ) is a diagonal matrix. This 
model has been discussed by lKuan and Led (120041 ) for testing a martingale difference hypothesis. 



Our proposed independence test can be applied to this nonlinear MA model, and the power 
values listed in Table [8] show that the proposed test performs well for this model. 

This discussion also indicates that the limit of the empirical spectral distribution of the 
nonlinear MA model (|5.8p is not the M-P law since the proposed test statistic is established on 
the characteristic function of the M-P law. 



5.5.2 Multiple ARCH(l) model 

Consider a multiple autoregressive conditional heteroscedastic (ARCH(l)) model of the form: 

W tj = Z tj y/ao + aiW?_ ld , t = l,...,p; j = l,...,n; (5.9) 

where zt = (Zn, . . . , Zt n ) is an n-dimensional random vector with i.i.d. elements, each of which 
has zero mean and unit variance, and u t = (Wti, ■ ■ ■ , Wt n )- 

For each j = 1, . . . ,n, ARCH(l) model (Wij,W2j, ■ ■ ■ ,W p j) is a martingale difference se- 
quence. ARCH(l) model has many applications in financial analysis. There exists no theoret- 
ical results stating that the LSD of the sample covariance matrix for ARCH(l) model is M-P 
Law, but from Figure [lj we can see that the LSD of the sample covariance matrix for ARCH(l) 
model is indeed M-P Law. A rigorous study is under investigation. For the ARCH(l) model, 
the proposed test can not capture the dependence of (u>i,u;2, • • • ,^ P ) directly, but we can test 
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the dependence of (u^ , w|, • • • , Since this test can tell us that (uf, u>|, ■ ■ ■ , <^>p) are not 
independent, naturally it can be concluded that (u>i, U2, ■ ■ ■ , are not independent either. 
Here, we take ceo = 0.9 and ol\ = 0.1. Table [9] shows the power values of our test for testing 
dependence of ARCH(l) model. 



5.5.3 Vandermonde matrix 

Consider the n x p vandermonde matrix V of the form 

/ 



V 



1 



-ILU2 



V 



-i(n— l)uii 



-i(n—l)u2 



~i(n—l)uj p 



(5.10) 



where cjj for i = 1, . . . ,p are called phased distributions and are assumed to be i.i.d on [0, 2tt). 
Then, the entries of V lie on the unit circle. Obviously, all the entries of the rows of V are not 
independent while the columns are independent. Denote the sample covariance matrix of V by 
D = V^V. 

Vandermonde matrices play an important role in signal processing and wi reless applications 



Ryan and Debban 



such a s direction of arrival estimation, pre-coding or sparse sampling theory. 
(J2009J) have established that as both n,p go to oo with their ratio being a positive constant, 
the limiting spectral distribution of D = V^V is not the M-P law. This result reminds us 
that the proposed test should be applied to capture the dependence structure of the rows of the 



matrix V. It is easy to see that, for any k = 1, 



1 and j = 1, . . . ,p, E(e 



-iku)j\2 



and 



E\e 



-ikujj 1 4 



1. The empirical power values given in Table [TU1 show that the proposed test works 



well in detecting dependence of Vandermonde matrices. 



6 Empirical analysis of financial data 

As an application of the proposed independence test, we test whether there is any cross-sectional 
independence among the daily closed stock prices of some relevant companies from New York 
Stock Exchange (NYSE) during the period 1990.1.1 — 2002.1.1. The data set is obtained from 
Wharton Research Data Services (WRDS) database. It consists of daily closed stock prices of 
632 companies at 5001 days. 

We randomly choose p companies from the total 632 companies in the data set. For 
each company i = 1,2, ...,p, there exist 5001 daily closed stock prices and denote it as 
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Si = (sii, Si2, ■ ■ ■ , Sj5ooi)- From the price vector Sj, we construct a time series of closed stock prices 
with length n: Xj = (sn, Si,i+50, Si,i+2x50, Sj,i+50n)- The construction of Xj,i = is 
based on the idea that the prices tend to be independent as the length of the time between them 
is large. Hence for each company i, the elements of stock price time series Xj are independent. 

The proposed test M n is applied to test cross-sectional independence of xi, X2, . . . , x p . For 
each (p,n), we randomly choose p companies from the database, construct the corresponding 
vectors Xi,X2, . . . , x p , and then calculate the P-value of the proposed test. Repeat this procedure 
100 times and plot the P-value graph to see whether the cross-sectional 'dependence' feature is 
popular in NYSE. 

From Figure [2l we can see that, as the number of companies p increases, more experiments are 
rejected in terms of the P-values below 0.05. When p = 30, n = 35, all the 100 experiments are 
rejected. This phenomenon is reasonable since as p increases, the opportunity that we choose 
dependent companies increases as the number of total companies is invariant. It shows that 
cross-sectional dependence exists and is popular in NYSE. This suggests that the assumption 
that cross-sectional independence in such empirical studies may not be appropriate. 

As comparison, we select 40 companies from the transportation section of NYSE and inves- 
tigate the dependence of daily stock prices of these companies from the same industry section- 
transportation section. We also perform the same experiments as above for those total 40 
companies. The P-value graphs of n = 10, p = 5, n = 15, p = 10 and n = 20, p = 15 are 
provided in Figure El Compared with the corresponding P-value graphs n = 10, p = 5, n = 15, 
p = 10 and n = 20, p = 15 in Figure [21 the daily closed stock prices of companies from the same 
section are more likely to be dependent. 



7 Conclusions 



This paper provides an approach for testing independence among a large number of high dimen- 
sional random vectors based on the characteristic function of the empirical spectral distribution 
of the sample covariance matrix of the random vectors. This test can capture various kinds 
of dependent structures, such as. MA(1), AR(1) model, nonlinear MA(1) model, ARCH(l) 



model and the general panel data mod e 



established in the simulation section. The conventional 



method (LRT proposed by lAndersonl (j 19841 )) utilises the correlated relationship between the 
random vectors to capture their dependence. This idea is only efficient for normally distributed 
data. It may be an inappropriate tool for non-Gaussian distributed data, such as martingale 
difference sequences (such as, ARCH(l) model), nonlinear MA(1) model, and the Vandermonde 
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matrix, which possess dependent but uncorrelated structures. The proposed test is not restricted 
to normally distributed data. Moreover, it is designed for a large number of high dimensional 
random vectors. An empirical application to closed stock prices of several companies from New 
York Stock Exchange highlights this approach. 



8 Appendix 



A.l Some useful lemmas 



Lemma 1 (Theorem 8.1 of lBillingsleyl (|1999l )). Let P n and P be probability measures on (C, (p). 
If the finite dimensional distributions of P n converge weakly to those of P, and if {P n } is tight, 
then P n =>• P. 



Billingsley 



(|1999l )). The sequence {X n } is tight if it satisfies these 



Lemma 2 (Theorem 12.3 of 
two conditions 

(I) The sequence {X n (0)} is tight. 

(II) There exists constants 7 > 0, a > 1, and a nondecreasing, continuous function F on 
[0, 1] such that 

E{\X n (t 2 ) - X n (hW} < \F(t 2 ) - F(h)\ a (8.1) 

holds for all t\,t2, and n. 

Lemma 3 (Continuous Theorem). Let X n and X be random elements defined on a metric space 
S. Suppose g : S — > S has a set of discontinuous points D g such that P(X E D g ) = 0. Then 



X n Ax^ g(X n ) A g{X). 



1.2) 



Lemma 4 (Complex mean value theorem (see Lemma 2.4 of lGuo and Nicolud (|2006l ))). Let f2 

be an open convex set in C. If f : f2 — >■ C is an analytic function and a, b are distinct points in 
Q, then there exist points u, v on L(a, b) such that 

Rei M^m ) = R<f{u% M /WlM) = lm(f\v)), (8.3) 
a — b a — b 

where Re(z) and Im(z) are the real and imaginary parts of z respectively; and L(a,b) = {a + 

t(b — a) : t £ (0,1)}. 



Lemma 5 (Theorem A. 44 of iBai and Silversteinl (|2009l )). Let A and B be two p x n complex 
matrices. Then, 



\F 



AA* 77BB 



1 



i 7 II < -rank(A-B) 
P 



(8.4) 
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A. 2 Proofs of main theorems 



Bai and Silverstein 



(|2004l ) established the remarkable central limit theorem for functional 



of eig envalues of A n under the additional assumption that -ElXnl 4 = 3 while 



Pan and Zhou 



(|2008l ) provid ed a supplement to th is theorem by eliminating the condition to some extent. From 



Theorem 4 of 



Pan and Zhoul (|2008l ) and (2.11) we can directly obtain the following lemma. 



Lemma 6. Under Assumptions 1 and 2, we have, for any positive integer k, 

cos(tix)dG n (x), . . . , J cos(tkx)dG n (x), J sin{t\x)dG n (x), . . . , J sin(tkx)dG n (x) S j (8.5) 
converges in distribution to the Gaussian vector of the form 

(V{h) + ^(ii), . . . , V(t k ) + tfi(tfc), Zl + S 2 (h), Z(t k ) + 6 2 {t k )), 
where 6i(t), 5 2 (t) are, respectively, defined as 

5i(t)=lim cos(tx)dp(F Cn ' Hn {x) - F Cn {x)), (8.6) 

n— >oo J 

S 2 (t) = lim sin(tx)dp(F Cn ' Hn (x) - F Cn {x)). (8.7) 

n— >oo J 

The means and covariances ofV(tj) and Z(tj) are specified in Theorem 1. 
Remark 3. When T n = I, the mean and varianc e of the asymptotic G auss ian distribution for 



Bai and Silverstein 



power functions f(x) = i r , Vr 6 Z + is calculated in \Pan and Zhou {(2008) and 
1 2004) - Hence the corresponding means and covariances for fi(x) = sintx and f 2 (x) = costx 
can be derived by Taylor series of sintx and costx. 

Proof of Theorem 1. Let t belong to a closed interval I = \T\,T 2 \. To finish Theorem 1, in 

view of Lemma [U and Lemma [6l it suffices to prove the tightness of { (</> n (t), ip n {t)) 

Thus it suffices to prove the ti ghtness of p(s n (t) — s(t)) . Repeating the same truncation and 



centralization steps as those in 



Bai and Silverstein! (|2004l ) . we may assume that 



\X{j\ < 5 n y/n, EXij — 0, E\Xij\ 2 



1, E\Xij\ A < oo. 



Set M n (z) = n[m F A n (z) — m F cn,H n (z)]. By the Cauchy theorem 

i r m 



2iri J z — x 



-dz, 



we have, with probability one, for all n large, 



e Ux dp(F A "(x) -F Cn ' Hn (x)) = d> e ltz M n {z)dz 



?Cn,H n 



2ni 



(8.9) 



(8.10) 
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The contour C involved in the above integral is specified as follows. Let 

C u = {x + iv : x £ [xi,x r ]}, 



.11) 



where vo > 0, x r is any number greater than limsup A max (T n )(l + \/c) 2 , %i is any negative 

n 

number if c > 1 and otherwise choose xi £ (0, limsup A m j n (T n )(l — \fc) 2 ). Then the contour C 

n 

is defined by the union of C+ and its symmetric part C_ with respect to the x-axis, where 



C, 



{xi + iv : v G [0, vq]} UC u U {x r + to : v G [0, uq]}- 



.12) 



From Lemma 1 and the argument regarding equivalence in probability of M n (z) and its 



truncation version given in Page 563 in 



Bai and Silversteinl (120041 ) and Lemma 3 we have 

D 

IC JC 
where M(z) is a Gaussian process, the limit of M n (z). 



\M n {z)\\dz\ 



\M{z)\\dz\, 



5.13) 



We conclude from Lemma H] that for any 5 > 

it\z it 2 z 



sup 

|*i— *a|<ff,*i,*aeJ ' 



< sup 

*a]<<5,*i,*ae/ 



( e «i* _ e «a*)M n (z)dz\ 



\J ' (Re(ize u ^)) 2 + (lm(ize u * z )) 2 5\M n (z)\\dz\ 



< K5 



\M n {z)\\dz\ 



D 



K5 



\M{z)\\dz\ 



as n — > oo, 



J.14) 



where t% and £4 lies in the interval [Ti,T2]> the last inequality uses (|8.13p and the fact that 
Re(ize ttsZ ), Im(ize zt4Z ) are bounded on the contour C and K (and in the sequel) is a constant 
number which may be different from line to line. 
By (|8.14|) . we have for any e > 0, 



P( sup 

v |*i-* 2 |«5,*i,* 2 e[o,i] 



(e itiz - e it2Z )M n (z)dz 



> e) < P(K5 



\M n (z)\\dz\ 



and 



lim lim sup P ( K5 

<5— >0 n-»oo 



\M n (z)\\dz\ 



> £ 



lim P[K6 

<5->0 



\M(z)\\dz\ 



Hence ([8T5]) and (f8TT6|) imply that 



lim lim sup P 



8^0 



|*i -ta 



sup 

<8,t!,t 2 £l 



it\z _ J,tiz 



)M n (z)dz 



> e 



> e 



0. 



> e 



0. 



$.15) 



.16) 



$.17) 



By Theorem 7.3 of lBillingslevI (|1999l ). J e ttx dp(F An (x) — F Cn ' Hn {x)) is tight. Moreover from 



the assumption we see that f e ltx dp(F Cn,Hn (x) — F Cn (x)) is tight by Lemma 4. 



□ 



23 



To prove Theorem 2, we first need to establish Lemma [7] below. To this end, write 

H n (x) = p[F s "(x) - F c "(x)\, (8.18) 

where S n is defined in (4.4). 

Lemma 7. Under the assumptions of Theorem 2, we have for any positive integer k, 

cos{t\x)dH n {x), . . . , J cos{tkx)dH n {x), J sin{t\x)dH n {x), . . . , J sin(tkx)dH n (x)^ (8.19) 

converges in distribution to a Gaussian vector (W(t\), . . . , W(tk), Q(ti), ■ ■ ■ , Q(tk)) whose mean 
and covariance function are given in Theorem 2. 

Proof of Lemma\7\ Using th e sam e truncation and centralization steps as those used in the 
paper by lBai and Silversteinl (|2004l ) . we may assume that 



Eij\ < Sn^/n, Esij — 0, E\si 



1, E\sij\ < oo. 



For the panel data model (4.1) proposed in Section 4, let 



The model can be written in the vector form as 



Vj = (vij,...,v pj Y , £j = (eij,...,e pj ) T , u= (— iti, . . . , -^pu p ) T , j = 1, 

VP 



1 



,n. 



.20) 



.21) 



Vj =£j+u, j = 1, ...,n. 

n 

We then define the sample covariance matrix by S n = ^ ^ v i v J • Moreover write 

3=1 



^ n 1 n 

V= — > V,-, £ = -V £j 



3=1 



3=1 



and 



n 1 n 1 n 

T>n = ~Yl £ i £ J> Sn = n ^ (Vj " ^ )(Vj " ^ )T ' Vn = n ~ 



£j - £ 



3=1 3=1 3=1 

Note that S n = T> n . The sample covariance matrix S n can be then expressed as 

S n = S n + vv T = T> n + vv T . 
By the conditions of Theorem 2 and the Burkholder inequality we have 



i=i 



■/? J 



7T 



.22) 



$.23) 



.24) 



.25) 
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which, together with Borel-Cantelli's Lemma, implies that 

X - as. 



u e 



0. 



Also, the conditions of Theorem 2 imply that 



u T u — > u. 



.26) 



Therefore, by the conditions of Theorem 2 and Theorem 2 of lPan and Zhoul (]201lh . we have, 
as n — > oo, 



A max (vv T ) = v T v = s T e + u T u + 2u T e ^4' c + u, as n — > oo. 
Furthermore, 



Jiang (|2004l ) proves that 

Kwx^Pn) ->■ (1 + \/c) 2 , a.S. US Tl Y OO 



and 



Xiao and Zhou 



(|2010l ) implies that, when c < 1 

A m i B (P„) — )■ (1 — \/c) 2 ) a - s - as ^ — )• oo. 



S.27) 



$.28) 



.29) 



By (|8.27p (|8.28p and (|8.29p . the maximal and minimal eigenvalues of S n satisfy with proba- 
bility one 



and 



limsupA max (S n ) < c + u+ (1 + \fc) 2 , 



liminf A min (S n ) > (1 - y/c)' 



.30) 



.31) 



As in the proof of Theorem 2, we obtain from Cauchy's formula, with probability one, for n 
large, 

P f I fi z ) 



p / f(x)d[F* n (x) -F c "(x)} 



2m 
V 

2-Ki 



z — X 



-dzd[F* n (x) - F Cn (x)} 



7 

f(z)dz 



1 



z — X 



d[F s "(x) - F Cn (x)] 



~ j f(z)(tr(S n - z^y 1 - pm Cn (z))dz, (8.32) 



where m Cn (z) is obtained from m(z) with c replaced by c n . The contour 7 is specified as follows: 
Let vq > be arbitrary and set 7^ = {/i + ivo, \i G [/x^, //,•]}, where \x T > c + u + (1 + -v/c) 2 and 
< (J,{. < 7( 0) i)(c)(l — \/c) 2 or fi£ is any negative number if c > 1. Then define 



7 + = + iv : v £ [0, v ]} U 7^ U + i« : v £ [0, i> ]} 



$.33) 
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and let 7 be the symmetric part of 7 + about the real axis. Then set 7 = j + U 7 . 



Set 



-1 
•-pi 1 



Then we have 



S-\z) = (S n - zl p r\ S~\z) = (S n - zl v 
V-\z) = (V n - zip)" 1 , D" 1 ^) = (D n - zip)- 1 . 

v T V~\z)v 



trSn 1 ^) ~ P m c n (z) = {trV n l {z) - pm Cn (z)) - 

l+V 1 V n (Z)V 



where we have used the identity 



(C + rr 



T\-l _ 



C^n^C" 1 



(8.34) 



(8.35) 



(8.36) 



l + r T C -i r ' 

where C and (C + rr T ) are both in vertible; and r 6 W. The first term on the right hand of 
(I8.35P was investigated in iPanl (|201ll ) . In what follows we consider the second term on the right 
hand of (15^51) . 

One may verify that 



(C + grv 



T\-l 



(8.37) 



l + gv r C-V 

where C and (C + qrv T ) are both invertible, q is a scalar and r, v £ MP. This, together with 
<K22\i and (fOBl . yields 



v T V-\z)v = e T V- 1 (z)s + 2u T V; i 1 (z)e + u T V; i 1 (z)u 



+ 2 \ ' + u t V- 1 {z)vl 



and 



v T P" 2 (z)v = e T V; 2 {z)e + 2u T V~ 2 (z)s + u T P" 2 (z)u 

e T Dn 2 (^ 2u^D- 2 (z) £ - ^ 2u r D- 1 (z) £ -^D- 2 (z)s ^ 



3.38) 



(8.39) 



(1 - £D n "^)e) 2 ^ 1 - e^D^e + (1 - £^D^(z)e) 2 
It is proved in Section 2.5 and (4.3) of IPanl (|201ll ) that as n — > 00, 

cm 2 (z) 



+ u T V- 2 {z)u. 



sup 

2:67 



e T B- 2 (z)s 



(1 + m{z)) 2 — cm 2 (z] 



0: 



(8.40) 



sup 

2G7 



I. p. 



0; 



and 



sup 

2S7 



i. p. 



(8.41) 



(8.42) 
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where we have also used an argument s imila r to (2.28) of iPanl (12011 )). 



By Lemma 2 of lBai. Miao and Pan! <j2007T ). (3.4) and (4.3) in[Pan 

as n — > oo, 



(2011 



sup 

ZG7 



u V n l (z)u — um(z] 



i.p. 



0. 



The next aim is to prove that 



sup 

ZG7 



u T T>„ 2 (z)u — u 



;dF MP (X) 



(\- z y 



i.p. 



and that 



sup 

zS7 



u T B~\z)e 



i.p 



0. 



)), and (4.5) we have, 
(8.43) 

(8.44) 

(8.45) 



Consider (|8.44p first. By the formula (18.36p . we have an expansion 



vi 1 V-\z)u = vi 1 D ? ; 2 (z)u + 



l-^D^i 



+ 



z)e 



(l-^D^(z)e) 2 



For any given z G 7, we conclude from Theorem 1 of IPanl (|201ll ) and Helly-Bray's theorem 



that 



u T T> n 2 (z)u — u 



MPtw l -P\ 



(\- z y 



;dF (X) 



as n — > 00. 



$.46) 



By the expansion of u T V n 2 (z)u and (j8.40j) - (|8.42p . to prove (|8.44|) . it suffices to prove the 
tightness of \k£\z) = u r D~ 2 (z)u -uf ^l^dF MP (A), z G 7} and |u- r D~ 2 (z)e, z G 7}. 



To this end, as in 



Bai and Silversteinl (120041 ). below introduce the truncated version of 



u T D n 2 (z)u. Define 7 r = {p r + iv : v G [n 1 p n ,Vo]}, 



{pi + iv : v G [n 1 p n ,v ]}, pi>0, 
{fit + iv : v G [0,v ]}, m < 0, 



$.47) 



where 



p n I 0, p n > n 6 , for some 6 G (0, 1) 



.48) 



Let 7+ = 71 U 7^ U 7 r and 7 n denote the symmetric part of 7+ with respect to the real axis. 
We then define the truncated process u T D^ 2 (z)u of the process u T D~ 2 (z)u for z = a + iv by 



U T Bn 2 (z)u 



u T D" 2 (z)u 



z G 7 n = 7„ U 7 n , 



^/D n - 2 (, n )u+^D„- 2 (, r2 )u M = p r ,v G [-n"V n ,n"V n ], 



.49) 



2p„ 



^u T D" 2 (z £l )u + £^u T D- 2 (^ 2 )u = w>0,»€ [-n~ Vn, n~Vn], 
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where z ri = p r + in p n , z r2 = \i r — in p n , = n e + in p n and ze 2 = pe — in 1 p n . We then 
have 

1 1 



sup 

zS7 



u T, D^ 2 (z)u-u i D; 2 (z)u <Kp n \\u\\\ 



+ ■ 



i. p. 



Amax(Dn) ~ Mr I I -^min(D n ) — [i[ \ 



0. (8.50) 



It is proved in Section 3 of iBai and Silversteinl (|2004l ) that, for any positive integer k and 

z 6 In U 7n> 

max(£'||D~ 1 (z)|| fc ) < ET. (8.51) 
It follows from independence between u and Sj,j = ,n that 

E\n T T>-\z)u- uJ J^dF MP (X)\ 



< E\u T V- 2 (z)u\ + \uj ^±-^ d F MP (X) 

< E\\u T \\ 2 E\\T>- 2 (z)\\ 2 + K<K, 



(8.52) 



which ensures Condition (1) of Lemma [2j Similarly, we can derive E\u T ~D- 2 (z)i\ 2 < K. 
Next, we prove condition (2) of Lemma [21 i.e. 



sup 



E\K^\ Zl )-K^{z 2 )\ 2 
\zi - z 2 \ 2 



< OO, 2=1,2. 



Note that 



i_ B i 



A (B — A)B , 



where A and B are any two nonsingular matrices. We then conclude that 

B~ 2 ( Zl ) - D^ 2 (z 2 ) = [z x - z 2 )n- 2 ( Zl )B-\z 2 ) + ( Zl - z 2 )D^ 1 (z 1 )D- 2 (z 2 ). 
Then 

k£\ Zi )-k£\z 2 ) 



(8.53) 



(8.54) 



(8.55) 



z\ - z 2 



u r D^ 2 (^ 1 )D- 1 (^)u + u T D^ 1 (^i)D; 2 (^)u 



{X - Zl) + {X - Z2) -dF MP {\). 



(\- Zl y(\- Z2 y 

As in (I8.52D . we can obtain 

.BI^D^^OD; 1 ^)^ 2 < K, B|u T D- 1 (z 1 )D; 2 (z 2 )u| 2 < K. 



(8.56) 



(8.57) 



Since /(A) = ^ s a continuous function when z\, z 2 E 7, the integral J f(X)dF MP (A) 



is bounded. This, together with (|8.57p . implies 

E\K n l \ Zl ) - K n l \z 2 )\ 2 



sup 

n, 21,2267 



< OO. 



\Zl - 22 1 
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(8.58) 



By (3.17) in iPan and Zhoul (l201ll ) and an argument similar to (I8.55p -( j8.58p we may verify 
that u T T)~ 2 (z)i is tight for z G 7. Summarizing the above we obtain (I8.44p . 

Consider ()8.45p now. From the last paragraph we see that it is enough to consider the 
pointwise convergence of u T D" 2 (z)e. As in (|8^9l) and (l830|) we may define the truncated 
process u T T>n 2 (z)e of the proces s u Y)„ 2 (z)e and the n prove that their difference tends to 
zero in probability. As in (4.3) of IPan and Zhoul (|201ll ) one may prove that for given z £ 7^ 
u T D r 7 2 (z)e ^ 0. 

From (j8738]) to (^4"5|) we have 



sup 

zS7 



<r> 1 . . , 1 + zm(z) , . \ 



and 



sup 

267 



v r p- 2 (,)v 



+ u 



■z 2 [(l + m(z)) 2 — cm 2 (z)} 
We then conclude from Slutsky's theorem that 



(A 



dF {Xj) 



t.p. 



$.59) 



.60) 



sup 

zS7 



V i P n 2 (z)v 4(l+m(2)) 2 -qn 2 (2) 



(A- 2 ) 



-dF (A) 



l + v r 2?r(z)v 



?jzm 2 (z) — 1 



.61) 



The arguments of Theorem 1 of lPanl ( 201ll ) show that the truncation version of (trP~ 1 (z) — 
pm Cn (2)) converges in distribution to a two-dimensional Gaussian process and that the dif- 
ference between (trT>~ 1 (z) — ym r ,„( z )) an d its truncati on ve r sion g oes to zero in probability 

)). Theorem 3 then follows 
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(see Page 563 of iBai and Silversteinl (|2004h and (2.28) of 
from ([HUH. (18.321) and (18.351) . Slutsky's theorem and Lemma 3 (one may refer to Page 563 of 
Bai and Silversteinl J2OOJ)). □ 



Proof of Theorem 2. As in the proof of Theorem 1, in view of Lemma [7] it suffices to prove the 
tightness of {p(s n (t) — s(t)) : t € /}. As before, write 



p(s n (t) - s(t)) 



p J e itx d[F Sn (x) - F Cn {x)} 



.62) 



where the contour 7 is specified in Lemma [71 
From the formula (|8.36p . we have 

trV- 1 (z)=trB- 1 (z) + 



$.63) 
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This, together with (I8.35|) . yields 

trS'Hz) -pm Cn (z) = trn- 1 (z)-pm Cn (z) + " Aj g _ n \\ (8.64) 

1 - e 1 D„ (z)e 1 + v 1 V n (z)v 

By (|8.32p and noting that M n (z) = trD~ 1 (z) —pm Cn (z), it is sufficient to prove the tightness 
of the following three terms: 

g n i(t) = -^-^e itz M n (z)dz, (8.65) 
9n2{t) = ~J e te l ^'ffij fe, (8.66) 

= "^7 / e te f3n- ^' (8 ' 67) 

The tightness of {<?ni(i) ■ t £ I = [T\, T2]} has been proved in Theorem 1. Next, via the same 
method adopted by Theorem 1, we prove the tightness of {g n i(t) : t £ I = [T\,T2]}, i = 2,3 as 
follows. 

By (pSlff|) . ([8H|) and Slutsky's theorem, we have 

i T ~D~ 2 (z)i cm(z) 



sup 

ZG7 



+ 



l-£ T D^(z)e z((l + m(z)) 2 -cm 2 ( 2 )) 
We conclude from (|8,68p . (|8.6ip and Lemma 3 that, as n — > 00, 
e T T>~ 2 {z)e 1 a .s. /1 cm(z) 



and 



^> 0. (8.68) 



^ f 1 77— ^ ^Hcfal (8.69) 

1 - e 1 D re (z)e Jj z[{l + m(z)) 2 - cm 2 (z)) 



V T P" 2 (Z)V , . a.s 



cm(z) 



l + v^-^v 1 ' 1 J 7 uzm 2 (z)-l 11 V ; 

By dSS2J, <K7Uh and the same proof as (gngj) to (|ST7|) . the tightness of {>„*(<) : t £ I = 
[Ti,T 2 ]}, i = 2, 3 can be derived. □ 
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Figure 2: P-values of the proposed test for daily closed stock prices from NYSE 
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*These are P-values for p companies from NYSE, each of which has n closed stock prices during the period 
1990.1.1 — 2002.1.1. The number of repeated experiments are 100. All the closed stock prices are from WRDS 

database. 
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Figure 3: P-values of the proposed test for daily closed stock prices from transportation sections 
of NYSE 



n=10, p=5 
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*These are P-values for p companies from transportation sections of NYSE, each of which has n closed stock 
prices during the period 1990.1.1 — 2002.1.1. The number of repeated experiments are 100. All the closed stock 

prices are from WRDS database. 
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Table 1: Empirical sizes and powers of the proposed test at significant level 0.05 for standardized 
normal distributed random vectors. 
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120 


0. 


056 





00(1 





091 


(1 


.059 


(1 


.002 


0.051 


0.045 


(1 


(175 


(1 


(178 





.078 


0.063 


30 





051 





062 





084 


(1 


.008 





.051 


0.081 


0.092 


(1 


.087 





115 





128 


0.148 


40 


0. 


041 





03(1 





082 


(1 


(155 


(1 


126 


0.183 


0.144 


(1 


208 


(1 


2(1!) 





27(1 


0.254 


50 





008 





081 





084 


(1 


.00(1 





151 


0.262 


0.329 


(1 


.420 





474 





4.32 


0.574 


60 


0. 


046 





027 


(1 


088 


(1 


.008 


(1 


154 


0.310 


0.455 


(1 


502 


(1 


896 





.733 


0.715 


70 





050 





044 


11 


018 


(1 


.110 





211 


0.379 


0.565 


(1 


747 





793 





84(1 


0.841 


80 





053 





04(, 





(131, 


(1 


102 





.233 


0.439 


0.687 


(1 


.892 





906 





944 


0.977 


90 


0. 


028 


0. 


(14(1 


(1 


054 


(1 


122 


(1 


22(1 


0.484 


0.752 


(1 


.878 


(1 


.963 





.980 


0.999 



*The powers are under the alternative hypothesis that the population covariance matrix is E = 0.057 p + 0.951 p l p . 
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Table 2: Empirical sizes and powers of the LRT at significant level 0.05 for standardized normal 
distributed random vectors. 



V 



n 




5 


10 


20 


30 


40 


50 


60 


70 


80 


90 


100 














Empirical sizes 










5 


0. 


005 


0.475 


0.495 

















(1 








10 





030 


0.115 


0.504 


0.454 























120 


0. 


047 


0.030 


0.522 


1 


0.500 




















30 





040 


0.032 


0.076 


0.873 


1 


1 

















40 





044 


0.034 


0.057 


0.148 


0.982 


1 


1 


1 











50 


0. 


057 


0.034 


0.044 


0.062 


0.299 


1 


1 


1 


1 


1 





60 





04!) 


0.053 


0.040 


0.065 


0.101 


0.578 


1 


1 


1 


1 


1 


70 


0. 


038 


0.060 


0.049 


0.054 


0.108 


0.229 


0.784 


1 


1 


1 


1 


80 





055 


0.040 


0.043 


0.072 


0.079 


0.147 


0.442 


0.950 


1 


1 


1 


90 


0. 


054 


0.061 


0.047 


0.033 


0.048 


0.129 


0.233 


0.691 


0.991 


1 


1 


100 


0. 


055 


0.048 


0.066 


0.069 


0.076 


0.085 


0.185 


0.402 


0.844 


0.997 


1 














Empirical powers 










5 


0. 


oor, 


0.504 


0.505 


























10 





033 


0.125 


0.506 


0.485 























120 


0. 


051 


0.049 


0.587 


1 


0.496 


(1 

















30 





056 


0.062 


0.121 


0.912 


1 


1 

















40 


0. 


061 


0.077 


0.139 


0.311 


0.990 


1 


1 


1 











50 


0. 


092 


0.090 


0.150 


0.206 


0.576 


1 


1 


1 


1 


1 





60 





088 


0.142 


0.220 


0.266 


0.420 


0.849 


1 


1 


1 


1 


1 


70 


0. 


1012 


0.170 


0.222 


0.320 


0.463 


0.691 


0.984 


1 


1 


1 


1 


SO 





133 


0.169 


0.246 


0.392 


0.478 


0.666 


0.883 


0.994 


1 


1 


1 


90 


0. 


1212 


0.190 


0.331 


0.404 


0.526 


0.666 


0.819 


0.983 


1 


1 


1 


100 





153 


0.232 


0.411 


0.550 


0.658 


0.713 


0.842 


0.952 


0.995 


1 


1 



*The powers are under the alternative hypothesis that the population covariance matrix is E = 0.05/p + 0.951-ply 
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Table 3: Empirical sizes and powers of the proposed test at 0.05 significance level for standardized 
gamma distributed random vectors. 



V 


11 




5 




10 


20 


30 


40 


50 


60 


70 


80 


90 


100 
















Empirical 


sizes 










5 





080 





078 


0.068 


0.059 


0.060 


0.063 


0.061 


0.047 


0.050 


0.061 


0.044 


10 





066 





075 


0.068 


0.054 


0.046 


0.051 


0.047 


0.048 


0.050 


0.042 


0.052 


12(1 





057 





058 


0.072 


0.062 


0.053 


0.054 


0.053 


0.052 


0.049 


0.051 


0.048 


30 





069 





069 


0.069 


0.076 


0.058 


0.056 


0.054 


0.047 


0.067 


0.066 


0.054 


40 





061 





04!) 


0.058 


0.040 


0.065 


0.048 


0.063 


0.065 


0.068 


0.047 


0.067 


50 





053 





054 


0.055 


0.057 


0.059 


0.048 


0.067 


0.066 


0.043 


0.059 


0.053 


60 





059 





052 


0.057 


0.060 


0.052 


0.067 


0.058 


0.064 


0.064 


0.061 


0.069 


70 





044 





050 


0.064 


0.055 


0.071 


0.054 


0.067 


0.064 


0.051 


0.077 


0.048 


80 





04.-, 





050 


0.061 


0.043 


0.071 


0.055 


0.071 


0.053 


0.056 


0.070 


0.060 


90 





041 





007 


0.034 


0.056 


0.049 


0.054 


0.050 


0.060 


0.047 


0.060 


0.058 


100 





070 





045 


0.059 


0.055 


0.047 


0.062 


0.069 


0.057 


0.056 


0.060 


0.061 
















Emp 


rical powers 










5 





334 





575 


0.853 


0.944 


0.983 


0.989 


0.998 


0.997 


0.999 


0.998 


0.998 


10 





513 





838 


0.979 


0.999 


1 


1 


1 


1 


1 


1 


1 


120 





721 





970 


0.999 


1 


1 


1 


1 


1 


1 


1 


1 


30 





834 





998 


1 


1 


1 


1 


1 


1 


1 


1 


1 


40 





014 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


50 





952 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


60 





991 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


70 





0012 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


80 





998 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


90 





999 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


100 




1 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 



*The powers are under the alternative hypothesis that the population covariance matrix is E = 0.05J P + 0.951 p l p . 



Table 4: Empirical powers of the proposed test at 0.05 significance level for MA(1) model. 



V 


n 




5 


10 


20 




30 


40 


50 


60 


70 


80 


90 


100 


5 





096 


0.097 


0.097 





111 


0.210 


0.214 


0.198 


0.221 


0.201 


0.223 


0.199 


10 





082 


0.090 


0.173 





227 


0.374 


0.715 


0.722 


0.809 


0.903 


0.913 


0.934 


20 





000 


0.165 


0.400 





507 


0.683 


0.822 


0.951 


0.994 


1 


1 


1 


30 





007 


0.121 


0.611 





733 


0.803 


0.986 


1 


1 


1 


1 


1 


40 





091 


0.321 


0.653 





938 


0.968 


1 


1 


1 


1 


1 


1 


50 





139 


0.416 


0.910 





956 


0.998 


1 


1 


1 


1 


1 


1 


60 





117 


0.412 


0.918 





004 


1 


1 


1 


1 


1 


1 


1 


70 





146 


0.535 


0.983 




1 


1 


1 


1 


1 


1 


1 


1 


80 





138 


0.620 


0.990 




1 


1 


1 


1 


1 


1 


1 


1 


90 





239 


0.805 


0.999 




1 


1 


1 


1 


1 


1 


1 


1 


100 





132 


0.761 


1 




1 


1 


1 


1 


1 


1 


1 


1 
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Table 5: Empirical powers of the proposed test at 0.05 significance level for AR(1) model. 



n 


5 


10 


20 


30 


40 


50 


60 


70 


80 


90 


100 


5 


0.130 


0.094 


0.264 


0.253 


0.181 


0.169 


0.191 


0.134 


0.144 


0.129 


0.113 


10 


0.167 


0.289 


0.482 


0.640 


0.668 


0.609 


0.563 


0.524 


0.427 


0.384 


0.321 


12(1 


0.230 


0.544 


0.878 


0.271 


0.954 


0.994 


0.999 


0.997 


0.996 


0.984 


0.974 


30 


0.205 


0.602 


0.993 


0.999 


0.671 


0.936 


1 


1 


1 


1 


1 


40 


0.344 


0.916 


0.998 


1 


1 


0.982 


1 


1 


1 


1 


1 


50 


0.541 


0.984 


1 


1 


1 


1 


1 


1 


1 


1 


1 


60 


0.571 


0.986 


1 


1 


1 


1 


1 


1 


1 


1 


1 


70 


0.672 


0.997 


1 


1 


1 


1 


1 


1 


1 


1 


1 


80 


0.673 


0.999 


1 


1 


1 


1 


1 


1 


1 


1 


1 


90 


0.882 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


100 


0.801 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 



Table 6: Empirical powers of the proposed test at 0.05 significance level for SMA(l) model. 



V 



n 


5 


10 


20 


30 


40 


50 


60 


70 


80 


90 


100 


5 


0.330 


0.431 


0.782 


0.999 


1 


1 


1 


1 


1 


1 


1 


111 


0.647 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2(1 


0.967 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


30 


0.962 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


40 


0.998 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


50 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


(ill 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


70 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


80 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


90 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


100 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 
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Table 7: Empirical sizes and powers of the proposed test at 0.05 significance level for the general 
panel data model. 



p 



IX 




■:> 




10 




20 




30 




40 




50 




60 




70 




80 




00 


100 




















Empir 


ical 


sizes 




















5 





038 





045 





(147 





()57 


o 


.058 





.067 





00!) 





(109 





072 





070 


0.074 


1(1 


0. 


041 





042 





(J40 





040 


n 


050 





050 





005 





047 





.068 





00.8 


0.069 


20 





035 





042 





(14!) 





.056 




052 





040 





003 





.043 





000 





055 


0.057 


30 





043 





04!) 





(143 





.059 





048 





068 





05!) 





057 





040 





055 


0.047 


40 


0. 


048 





052 





(143 





.060 


n 


057 





040 





049 





()54 





040 





058 


0.061 


50 


0. 


057 





048 





040 





058 


n 


052 





055 





048 





.049 





(150 





(140 


0.041 


60 


0. 


058 





056 





055 





048 


n 


047 





045 





.053 





(100 





058 





.049 


0.050 


70 





062 





060 





05!) 





056 


[) 


04!) 





057 





040 





008 





052 





036 


0.043 


80 


0. 


071 





003 





007 





047 





048 





058 





050 





()44 





05!) 





057 


0.055 


90 





005 





068 





065 





()48 





.053 





.048 





050 





048 





048 





000 


0.060 


100 





037 





046 





03(i 





.0.35 





045 





043 





045 





.050 





048 





047 


0.055 






















Empirical p 


owcrs 
















5 





150 





238 





.345 





.417 





.484 





52!) 





540 





.615 





Oil 





.668 


0.692 


10 


0. 


125 





247 





452 





526 





568 





633 





669 





7.37 





705 





751 


0.800 


120 


0. 


206 





343 


(1 


493 





.615 





073 





752 





788 





.813 





.860 





.864 


0.876 


30 


0. 


111 





4(14 





535 





.684 





757 





756 





.855 





.875 





882 





.909 


0.953 


40 


0. 


308 





393 


(1 


005 





.698 





780 





820 





878 





.898 





044 





.959 


0.953 


50 





207 





450 





603 





718 





815 





.889 





923 





.9.88 





966 





973 


0.980 


60 


0. 


268 





430 


(1 


594 





780 





820 





018 





913 





.926 





974 





970 


0.984 


70 





144 





434 





04!) 





708 





.888 





.883 





.944 





.908 





971 





982 


0.996 


80 


0. 


171 





.454 





078 





.796 





872 





021 





.938 





.967 





98!) 





992 


0.995 


90 





204 





.431 





083 





.,834 





874 





916 





003 





.985 





988 





994 


0.994 


100 





291 





.398 





087 





.836 





.884 





0.31 





073 





.987 





992 





994 


1 



Table 8: Empirical powers of the proposed test at 0.05 significance level for nonlinear MA model. 



V 



11 




5 




10 




20 




30 




40 




50 




60 




70 




80 




90 


100 


5 


0. 


033 





023 





004 





.008 





.005 





.007 





008 





.007 





.00!) 





014 


0.070 


10 


0. 


018 





.543 





374 





220 





108 





.037 





Oil 




















20 





804 





744 





703 





.614 





581 





.511 





447 





.340 





300 





257 


0.216 


80 


0. 


854 





.815 





777 





.779 





780 





740 





002 





.002 





507 





579 


0.555 


40 





878 





842 





856 





.856 





845 





.825 





77!) 





770 





772 





702 


0.698 


50 





884 





80!) 





.868 





.884 





804 





.888 





800 





.875 





80!) 





.828 


0.820 


00 





892 





912 





.882 





.904 





020 





928, 





033 





.900 





892 





892 


0.882 


70 





890 





.934 





906 





.927 





034 





921 





952 





.925 





943 





917 


0.926 


80 


0. 


936 





880 





025 





921 





952 





950 





943 





943 





958 





.958 


0.936 


90 





922 





925 





03!) 





.9.85 





958 





959 





955 





.982 





962 





957 


0.954 


100 


0. 


926 





.935 





.920 





.937 





952 





.964 





.970 





.975 





.978 





.970 


0.965 



40 



Table 9: Empirical powers of the proposed test at 0.05 significance level for ARCH(l) model. 



p 



11 




5 




10 




20 




30 




40 




50 




60 




70 




80 




90 


100 


5 


0. 


148 





093 





(173 





076 





058 





.099 





11.3 





(188 





106 





112 


0.141 


10 





454 





399 





315 





.216 





128 





121 





112 





121 





09.3 





095 


0.092 


20 





401 





428 


11 


508 





464 





444 





.468 





4.36 





.338 





356 





.337 


0.302 


30 


0. 


12712 





364 





616 





712 





753 





759 





79.3 





743 





.638 





.639 


0.598 


40 





2.32 





357 





572 





.823 





790 





874 





915 





.885 





876 





.855 


0.860 


50 


0. 


12212 





339 





022 





.757 





891 





.969 





9.17 





975 





967 





957 


0.975 


60 





216 





.448 





617 





.862 





901 





976 





983 





.987 





992 





997 


0.997 


70 





200 





339 





.592 





.864 





931 





982 





99.3 





.998 





997 





998 


0.996 


80 


0. 


104 





376 


(1 


.566 





824 





950 





.967 





997 




1 




1 





999 


1 


90 





184 





.456 





721 





.839 





960 





995 





999 




1 





999 




1 


1 


100 


0. 


128 


0. 


306 





802 





.859 





.934 





.992 




1 




1 




1 




1 


1 



Table 10: Empirical powers of the proposed test at 0.05 significance level for Vandermonde 
Matrix. 



11 




10 




20 




30 




40 




50 




60 




70 




80 




90 


100 


120 


10 





.180 





197 





192 





169 





202 





211 





190 





185 





171 





177 


0.240 


20 





.309 





332 





356 





327 





291 





808 





.301 





295 





.321 





24.8 


0.478 


30 





.824 





.433 





473 





.413 





461 





408 





445 





.395 





.368 





.897 


0.606 


40 





.458 





512 





527 





546 





.533 





490 





518 





498 





450 





457 


0.655 


50 





59.3 





437 


I) 


549 





571 





614 





569 





577 


0. 


566 





565 





537 


0.764 


60 





.504 





.538 





551 





567 





616 





662 





.588 





581 





572 





607 


0.744 


70 





548 





526 





560 





.627 





.668 





641 





.694 





707 





641 


0. 


678 


0.741 


80 





550 





545 





589 





.633 





712 





719 





69.3 





768 





729 





749 


0.805 


90 





.589 





544 


I) 


596 





.667 





695 





712 





74.3 





754 





738 





728 


0.807 


100 





.464 





549 





61(1 





.645 





704 





757 





772 





751 





752 





808 


0.928 


120 





.6.3.3 


() 


660 


(1 


736 





.737 





.759 





855 





854 


0. 


909 





.960 





999 


1 
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